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Abstract. Freeze out of particles across a space-time hypersurface is dis- 
cussed in kinetic models. The calculation of final momentum distribution of 
emitted particles is described for freeze out surfaces, with spacelike normals. 
The resulting non-equilibrium distribution does not resemble, the previously 
proposed, cut Jiittner distribution, and shows non-exponential p t -spectra sim- 
ilar to the ones observed in experiments. 
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1. Introduction 

Continuum and fluid dynamical models, due to their simplicity are very popular 
in heavy ion physics, because they connect directly collective macroscopic matter 
properties, like the Equation of State (EoS) or transport properties, to measurables. 
In contrast to kinetic, Monte Carlo models, in continuum models the evalua- 
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tion of measurables is a problem in itself. Particles, which leave the system and 
reach the detectors, can be taken into account via source (drain) terms in the 
4-dimensional space-time based on kinetic considerations, or in a more simplified 
way via freeze out (FO) or final break-up schemes, where the frozen out particles 
are formed on a 3-dimensional hypersurface in space-time. Such freeze out de- 
scriptions are important ingredients of evaluations of two-particle correlation data, 
transverse-, longitudinal-, radial-, and cylindrical- flow analyses, transverse momen- 
tum and transverse mass spectra and many other observables. 

The general theory of discontinuities in relativistic flow was not worked out 
for a long time, and the 1948 work of A. Taub jjj discussed discontinuities across 
propagating hypersurfaces only (which have a spacelike unit normal vector da^-, 
da^da^ = — 1). Events happening on a propagating, (2 dimensional) surface belong 
to this category. 

Another type of change in a continuum is an overall sudden change in a finite 
volume. This is represented by a hypersurface with a timelike normal (da^da^ = 
+1). If one applies Taub's formalism to freeze out surfaces with timelike normal 
vectors, one gets a usual Taub adiabat, but the equation of the Rayleigh line will 
yield imaginary values for the particle current across the front. Only in 1987 Taub's 
approach was generalized to both types of surfaces , making it possible to take into 
account conservation laws exactly across any surface of discontinuity in relativistic 
flow. This approach also eliminates the imaginary particle currents arising from the 
equation of the Rayleigh line. When the EoS is different on the two sides of the 
FO front these conservation laws yield changing temperature, density, flow velocity 
across the front. Based on conservation laws and on simple kinetic considerations 
the idealized surface freeze out models were completely worked out in refs. ||,f|,|6|,D ■ 



2. Conservation laws across idealized freeze out fronts 

The FO hypersurface is an idealization of a layer of finite thickness (of the order of 
a mean free path or collision time) where the frozen-out particles are formed and 
the interactions in the matter become negligible. The dynamics of this layer can be 
described in different kinetic models as Monte Carlo models or four-volume 

emission models [p^~|— p~5| . The zero thickness limit of such a layer is the idealized 
FO surface. Kinetic models for hadronic degrees of freedom indicate that such 
an idealization is meaningful only for collisions of massive heavy ions like Au+Au 
or Pb+Pb [To) . If we include quark-gluon plasma in our reaction model with 
rapid final hadronization which coincides with freeze out |l(||l7|], the applicability 
of idealized surface freeze out description becomes even better. 

In general the energy - momentum tensor changes discontinuously across this 
idealized hypersurface. If the flow is not orthogonal to this surface the four-vector of 



the flow velocity will also change across this surface p, 16 , lq| . The FO discontinuity 
has a timelike normal in most cases, and so the method for the description of timelike 
detonations and deflagrations j| should be used (see also @,|l8|-|l|). 
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If da^ is the normal vector of an element of the FO hypersurface (c£<7 M is a unit 
vector parallel to efcr M ), the observable triple differential cross section is obtained by 
integrating the FO current over the whole FO hypersurface. This is described by 
the Cooper-Frye formula 



= J fFo{x,p;T,n,u") p^da^ , (1) 

where fFo(xiP]T 1 n,u u ) is the post FO phase space distribution of frozen out par- 
ticles which is not known from the fluid dynamical model. When we apply this 
formula we have to ensure || that: (i) the particles are really leaving the surface 
outwards: 

f FO (x,p;T,n,u u ) = 0(p M da-^)f FO (x,p; T, n, u") , (2) 
and (ii) the conservation laws and entropy condition are satisfied (l], |) : 

[iV" da = , [T^ da^] = and [S" da M ] > , (3) 

where [A] = A — Aq, and the pre FO baryon and entropy currents and energy- 
momentum tensor are Nq, Sq and Tq V , while the post freeze out quantities are, 

and T'"', In numerical calculations the local freeze out surface can be determined 
most accurately via self-consistent iteration j|;[22| . This fixes the parameters of our 
post FO momentum distribution, fpo(x,p; T, n, u v ) if the shape of the distribution 
is known. 

A simple approach is to assume that fpo 1S a cu t Jiittner distribution ^,^| if we 
have a FO hypersurface with spacelike normal. It was shown that one can formally 
solve the freeze out problem with this ansatz . 



3. Freeze out distribution from kinetic theory 

At the same time in a kinetic gas model it was studied how we can obtain a non- 
equilibrium post FO distribution across a FO hypersurface with a spacelike normal 
. The result indicated that actually the cut Jiittner distribution is not a realistic 
choice. Apart of the non-physical shape by the sharp momentum cut, the simplified 
model which recovered the cut Jiittner shape led to other nonphysical consequences 
also: The freeze out was not complete in the model even for an asymptotically large 
thickness of the FO layer. Thus, the simple kinetic model || had two unsatisfactory 
features: (i) it did not achieve complete freeze out, and (ii) it resulted in exponen- 
tially weakening freeze out with a FO layer thickness which tends to infinity. A 
solution to the first problem has also been suggested. 

Here we present solutions to both problems by making additional improvements 
of the model used in ||. The two extended models are still very superficial and 
should not be considered as final and physically realistic solutions of the freeze 
out problem. For transparency we perform and demonstrate the modifications 
separately, although both of these could be included in a combined model. This 
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way we can see which improvement (consideration of which physical process) cures 
problems. 

Let us assume an infinitely long tube with its left half (x < 0) filled with nuclear 
matter and in the right vacuum is maintained. We can remove the dividing wall at 
t — 0, and then the matter will expand into the vacuum. By continuously removing 
particles at the right end of the tube and supplying particles on the left end, we can 
establish a stationary flow in the tube, where the particles will gradually freeze out 
in an exponential rarefaction wave propagating to the left in the matter. We can 
move with this rarefaction wave or front, so that we describe it from the reference 
frame of the front (RFF), where the rarefaction is stationary. 

In this frame, we have a stationary supply of equilibrated matter from the 
left, and a stationary rarefaction front on the right, x > 0. We can describe the 
freeze out kinetics on the r.h.s. of the tube assuming that we have two components 
of our momentum distribution []l2|-|l4]], ff ree (x,p) and fi n t{x,p). However, we only 
assume that at x = 0, ffree vanishes exactly and fi nt is an ideal Jiittner distribution 
(supplied by the inflow of equilibrated matter), while fi n t gradually disappears and 
ffree gradually builds up as x tends to infinity. 

Let us recall || the most simple kinetic model describing the evolution of such 
a system. Starting from a fully equilibrated Jiittner distribution the two compo- 
nents of the momentum distribution develop according to the coupled differential 
equations: 



COS u ~* 

d x f mt (x,p)dx = -Qip^da^) p f mt (x,p)dx, 



COS u ~* 

d x ff r ee(x,p)dx = +0(^0^) - P f mt (x,p)dx. (4) 

Here cos 6*^ = ^ in the RFF frame. The interacting component, fint, will 
deviate from the Jiittner shape and the solution will take the form: 



fint(x,f) = fjuttner(x = 0,p)exp 



-Q(p»dv»f-^x 



(5) 



This solution is depleted in the forward p-direction, particularly along the x-axis. 
Inserting it into the second differential equation above, leads to the freeze out solu- 
tion: 



ffree(x,p) = fjuttner(x = 0,p)<l- exp 



cos 9f; 



-e(p»da n )^^x 

,(,) 



fjuttner(x =0,p)~ fint{x,p) ■ 

At x — ► oo this distribution will tend to the cut Jiittner distribution introduced 
in the previous section. The remainder of the original Jiittner distribution survives 
as fint, even if x — > oo. In this model the particle density does not change with 
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x, barely particles moving faster than the freeze out front (i.e. p^da^ > 0) are 
transferred gradually from component fi n t to component // ree . This is a highly un- 
realistic model, indicating that rescattering and re-thermalization should be taken 
into account in fi nt . This would allow particle transfer from the " negative momen- 
tum part" (i.e. p^da^ < 0) of fi n t to // re e, which is not possible otherwise. 

4. Freeze out distribution with rescattering 

The assumption that the interacting part of the distribution remains the distorted 
(after some drain) Juttner distribution, is of course highly unrealistic. As suggested 
in ref. j^] rescattering within this component will lead to re-thermalization and 
re-equilibration of this component. Thus the evolution of the component, fi n t is 
determined by drain terms as much as the re-equilibration. Here, we repeat the 
introduction of the model with rescattering, and present it more generally for both 
positive and negative flow parameter velocities. This turned out necessary because 
during the FO process the interacting component will acquire negative velocity 
parameters, even with large positive initial velocities. We also introduce the m = 
limit of the suggested solution which enables us to present a transparent and not 
completely numerical solution. 

If we include the collision terms explicitly into the transport equations (Q) in 
general case leads to a combined set of integro-differential equations. We can, how- 
ever, take advantage of the relaxation time approximation to simplify the description 
of the dynamics. 

Then the two components of the momentum distribution develop according to 
the coupled differential equations: 

d x fi n t(x,p)dx = -<d(p tl da t _ l )^£f mt (x,p)dx+ 

(7) 

+ lfeq(x,p) - fint(x,p)] jrdx, 

d x ffree{x,p)dx = +@(p» da Al ) £2 ^ f fint(x , p)dx . (8) 

The interacting component of the momentum distribution, described by 
eq. (0), shows the tendency to approach an equilibrated distribution with a re- 
laxation length A'. Of course due to the energy, momentum and conserved particle 
drain, this distribution, f eq (x,p) is not the same as the initial Juttner distribution, 
but its parameters, n eq (x), T eq (x) and u% q (x), change as required by the conserva- 
tion laws. 

Conservation Laws. In this case the change of the conserved quantities caused 
by the particle transfer from component int to component free can be obtained in 
terms of the distribution functions as: 

dN? = ~ [ tZpPQbf&p) cos 8 p ~f mt (x,p) (9) 
A J Po 
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Fig. 1. The temperature and flow veloirty of the interacting component in the 
rest frame of the freeze out front (RFF) for baryonfree and massless gas, calculated 
according to the coupled set of equations ( f20|). We assume an initial temperature 
of T = 130 MeV to fit the situation in refs. 1^,0. 
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and 

dTr = _dx r fP^Q^M) cos%/ int (x,p). (10) 
A J Po 

If we do not have collision or relaxation terms in our transport equation then the 
conservation laws are trivially satisfied. If, however, collision or relaxation terms 
are present, these contribute to the change of T^ v and iV M , and this should be 
considered in the modified distribution function fmt{ x ,p)- 



4.1. Immediate re-thermalization limit 

We can get from the equation (fy the general solution for fi n t(x,p) in the RFF 
frame: 



+i£/ eg (*>>^)^W}. (11) 

As a first approximation to this solution let us assume that A' <C A, i.e. re- 
thermalization is much faster than particles freezing out, or much faster than pa- 
rameters, n eq (x), T eq (x) and u£ q (x) change. Eq. ( pi] ) then leads to 

fint(x,f) ^ f eg (x,p), for A'«A. (12) 

For f eq (x,p) we can assume the spherical Jiittner form at any x including both 
positive and negative momentum parts with parameters h(x), T(x) and u^ FG (x), 
where u^ FG (x) is the peak velocity of the Jiittner gas (which is the same as the 
flow velocity of the non-cut Jiittner gas), and h{x) is the proper density i.e. the 
density in the frame moving with u^ FG (x) 

In this case the change of conserved quantities due to particle drain or transfer 
can be evaluated for an infinitesimal dx. We assume that the 3-flow is normal to 
the freeze out surface. The changes of the conserved particle currents and energy- 
momentum tensor in the RFF, using the notation of ref. J(|, b eqs. (| 10) are given 

by 
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p x [100 MeV] 

Fig. 2. The freeze-out distribution, ff ree (x,p), in the rest frame of the freeze 
out front (RFF) calculated according to eqs. (po|). A, B and C correspond to 
x = 0.2A, 3A, 100A respectively and u^ FG \ x= o = (1,0,0,0). The numbers in the 
contours are in arbitrary units. The distribution is asymmetric and elongated in 
the freeze out direction, x. This may lead to a large-pt enhancement, compared 
to the usual Jiittner assumption used in most previous calculations as a freeze out 
distribution. Note that fi n t(x,p) does not tend to the cut Jiittner distribution in 
the limit x — > oo. 



dN° = -fj^ {aK x {b) + jb(3v 2 - 1) 7 2 [(1 + j)KM - M«, &)] 
+ n v 2 b 2 [(l + j)K (a) - JC (a, b)} + 2jv 3 j 3 (b + l)e~ b } 

=S -ff 7 (l +U )2, 

dN?= -f J -^{~2K (b) + 2jv 1 2 e- b lv 2 1 2 (b+l)-v 2 b-l} (13) 
+j(2 + v^ 2 b 2 )[(l + j)Ko(a) - /C (o, b)} + 
+jv 2 (3v 2 - l) 7 3 i[(l + j)ifi(fl) - JCi(a,b)]} 
^ [- 2 «7 2 (1 + u) - 2in(l - «) + u 2 7 4 (l + «) 2 ] , 
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and by the expressions 
dT? a = 



dx 
" A 



dTP* 



dT? v 



^ {^(6) +jv^e- b [(l+3v^A(b)- 

-(2 + v 2 b 2 )(b + 1) + v 4 (l + ^) 7 2 & 3 ] + a 2 K (b) 
+jv 2 1 2 b 2 (3 + v 2 ) [(1 + j)K 2 (a) - IC 2 (a, b)} 
+j(v 2 b 2 -v 2 - l)jb [(I +j)K 1 (a) - ld(a,b)] 

-jb 2 [(l+j)K (a)-IC (a,b)}} 
~^^[l-u 2 + 2u 1 \l + uf- 1 2 {l+u) 2 } , 

dx 0(1 + 3v 2 )v 2 ^b 2 [(1 + j)K 2 (a) - lC 2 (a, b)} 



A Ajv' 6 ~{' 6 



+jv 2 j 5 e h 



v(v 2 + 3)A(b) - a 2 bv 3 + (1 + 3i; 2 )^- - va 2 



+jv^b 3 [(l+m, (a) - /Ci (a, 6)] } - ^diV° 



-f ^ 7 2 (1 + M )3 , 



dZT = -f^i {jv 4 (3 + v 2 )j% 2 [(1 + j)K 2 (a) - JC 2 (a, b)} 
+jv e j 3 b 3 [(j + l)Ki(a) - /Ci(a,&)] +^ 3 7 6 e-" 



t>(3+n 2 )(H 



da; nT 



+ 



a 2 (« 4 fe-l) + (3t; 2 + l)A(fe)]}-|^ 



diV? 



(14) 



- w [2u 3 7 6 (l + «) 3 - 3« 2 7 4 (1 + «) 2 + 
+6u7 2 (l+u) + 6;n(l-u)] , 

-X W W^W + !) & [(1 +WM &)] 
-jt, 4 7 2 & 2 [(1 +i)Xo(a)-/C (a,6)] + 6« 2 ^i(6) 
-2^ 3 7 4 (&+l)e- b } + ;|^f 
^ -TW [2uV(l+ U ) 2 - 6 U7 2 (l+«) + u 2 - 2/n(l-«)] , 
and dT" = dT?". Note that in RFF the flow velocity of the re-thermalized 
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component is ' RFG {x) — j(x) (l,u(x),0, 0)\rff, where 7 = l/Vl — 
also use the notation v = \u\ and j — sgn(u). 

The new parameters of distribution fi n t, after moving to the right by dx can be 
obtained from dN^ and dT^ v . The conserved particle density of the re-thermalized 
spherical Jiittncr distribution after a step dx is: 



hi(x + dx) — hi{x) + dfii(x) = J N^(x + dx)Ni^(x + dx) , 

where the expressions are invariant scalars. The differential equation describing the 
change of the proper particle density is: 

dhi(x) = uf RFG (x) dN iifl (x) . (15) 

Although this covariant equation is valid in any frame, we can calculate it in the 
RFF, where the values of dN^ were given above in eq. (|l3|). 

For the re-thermalized interacting component Eckart's flow velocity is the ve- 
locity of the RFG, which changes with x, so we denote this frame as RFG(ir). 
The velocity of this frame decreases with decreasing x due to the particle drain at 
positive momenta. For the spherical Juttner distribution the Landau and Eckart 
flow velocities are the same, uf E RFG (x) — u% L RFG (x) — RFG (x). Thus we can 
evaluate the flow velocity v% RFG (x + dx): 



<rf G (i + dx) =N?(x + dx)/^N~, 
which leads to the following covariant expression 



where A^(x) = g^ v — RF q{x) u" rfg (x) is a projector to the plane orthogonal 
to it£ RF q{x). This equation is valid in any reference frame, nevertheless we know 
the four-vectors on the r.h.s. in the RFF explicitly. Then the new Eckart flow 
velocity of the matter is u? E RFG {x + dx) = u^ RFG (x) + du^ E RF q(x). 

To get the temperature and the change of Landau's flow velocity, we analyze 
the change of the energy momentum tensor. Before the particle drain the energy - 
momentum tensor at x in the RFG is diagonal, T^ u {x) — diag(ei, Pj, Pj, Pi)\ RF G{x) 
while in the RFF 1f{x) = [(e, + P t ) < flFG < ifiFG (x) -Pig""] \rff („) ■ Adding 
to this the drain terms, dT^ v {x), arising from the freeze out while we move to the 
right by dx, yields T- 11 \x + dx) which will not be diagonal in the RFG(x) and the 
pressure part will not be isotropic. We can Lorentz transform this to another frame 
which diagonalizes Tjf (x + dx). This means to find the Landau flow velocity of 
the new system, L RF q(x + dx) in the original RFG(ir). After a straightforward 
diagonalization, a somewhat tricky algebra and neglecting second and higher order 
terms we arrive at the covariant expression^ 

d K,L,RF G ( X ) = T p ■ ( 17 ) 

c-i 1 Lj 
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Fig. 3. The same as Fig. g, just u RFG \ x =o — 7(1,0.5,0,0). We can see that 
particles are freezing out faster if we start with positive initial flow velocity. 



Although, for the spherical Jiittncr distribution the Landau and Eckart flow ve- 
locities are the same, the change of this flow velocity calculated from the baryon 
current and from the energy current are different in general 

^ U i,E,B.FGi X ) 7^ dUi L RFG (x) . 

This is a clear consequence of the asymmetry caused by the freeze out process as it 
was discussed in ref. || . 

In the special case of the massless limit we can calculate du 1 ^ E RF q{x) and 



duf 



dul L RF g( x ) m the RFC In this frame = (1, 0, 0, 0), and according to ( pL 6| ) and 
( |l7| ) we can get: du® E — du® L — 0, du\ E — du\ L — 0, du\ 

dN'f dNf 



and 



dT 



xO 



dT, x0 



P T? 



Here the last values of both equations are in the RFC Thus we can see that 



.H) 



duf E — ^duf L , in the massless limit. Using the definitions of the dN x and dT 
from equations (p|),(|lC|) and N° and T 00 from section[]2., it is easy to see, that in 
the limit m — > 0, (and p° = p) the calculation leads to the integral: 



j. 4 „ dx 
du x iyE = -dul L = - — 



dcosO 



cos ft cost/ — u) 



ucos( 



8(tt — cos ft) , 
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P x [GeV] 

Fig. 4. The local transverse momentum (here p x ) distribution for baryon free, 
massless gas at p y — 0, x = 100A and To = 130 MeV. The transverse momentum 
spectrum is obviously curved due to the freeze out process, particularly for large 
initial flow velocities. The apparent slope parameter increases with increasing trans- 
verse momentum. This behavior agrees with observed pion transverse mass spectra 
at SPS @,|2| . 

where u is defined via u^ RFG = j(x)(l, u, 0, 0) (in RFF), and cos# is the polar 
angle of the emitted particle in the RFC (The polar emission angle in RFF, given 
earlier, can be expressed as cos 9p — (cos 6 — u)/(l — u cos 9) .) 

This problem does not occur for the freeze out of baryonfree plasma, and we 
have only duf* L . 

The last item is to determine the change of the temperature parameter of fi n t ■ 
From the relation e = u^T^Uy we readily obtain the expression for the change of 
energy density 

de l (x) = u^i >RFG (x) dT^ v (x) u v . MRFG (x) , (18) 

and from the relation between the energy density and the temperature (e.g. Chapter 
3 in ref. 0), we can obtain the new temperature at x + dx. Fixing these parameters 
we fully determined the spherical Jiittner approximation for fi nt . 

The application of this model to the baryonfree and massless gas gives the 
following coupled set of equations: 

dx 4ers b 5 

(19) 
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^ = --=-2— \t» v - U^U a T GV \ U v . 

Here we use the EoS, e = crg^T 4 , and the definition (ff Ml/ = —dxr^T 4 , where 
(iT^ are given by eqs. (M), and so x is measured in units of A. 

The results of numerical calculation are displayed in Figs. |l| - 1|. The velocity 
and temperature of the interacting component are gradually decreasing due to the 
loss of particles carrying most momentum and thermal energy. Thus, we see that 
T -> 0, when x — ► oo . So, fi n t {x,p) — ^■Kh) i ^pKa* — P" u v)/T\ — > 0, when x — > oo. 
Thus, all particles freeze out in the model with rescattering. 

Now we can find the distribution function for the noninteracting, frozen out part 
of particles according to equation (ph . The results are shown in Fig. || - Fig. ||. We 
would like to note that now ff ree (x,p) does not tend to the cut Jiittner distribution 
in the limit x — > oo, it has a smooth, anisotropic shape. Most importantly the slope 
in the FO direction is not exponential, resembling recent experimental data. 

5. Volume emission model 

In this section we demonstrate an improvement which yields complete FO in a layer 
of finite width. We calculate the kinetic freeze out distribution based on four volume 
emission models |I2-|l4j . In order to illustrate the physical mechanism of this freeze 
out process, let us study again a simple one-dimensional flow MM. We suppose an 
infinite tube where a stationary flow of a fluid is supplied from the left (a; < 0) , so 
that the freeze out occurs for the positive direction of x. 

In the four volume emission model, we introduce the basic quantity, the so- 
called escape probability 

V{r,t,p) = e ■>* , (20) 

where n is the total density in the calculational frame, v = p/E, the velocity of the 
particle, a, the total cross section and v re i, the relative velocity. 

We can understand this assumption in the following way: let us think for sim- 
plicity of an ideal gas of hard balls with radius R e ff, so that the collision cross 
section becomes (To = AitR^,^ . A particle is frozen out at time t if it will not collide 
any more starting from this time. The probability to have no collisions, or to freeze 
out, is described by a Poisson distribution 

P = W(0) = e~ pN , (21) 

where p is the probability of a collision, and N is the total number of particles. 
Then, in the hard ball approximation pN is the number of particles which our 
particle can meet on its way 



/oo 
a v re in(r+ p/Et,t)dt , 



(22) 
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where n(x) is the particle density in the calculational frame n(x) = N°(x). This 
leads us to the expression (^) , where the interaction can be included through hard 
ball cross section a = <To or the effective cross section tr = <r e / / • 

For a stationary one dimensional case, we can express the escape probability, 

V, as c 

r e -x><*^, cose>o 

V(r,t,p)^V(x,cos6) = i (23) 

( e -/," %n W^, cosfl<0 

where cos (9 = — . Let us rewrite this result in the following form 

V(x, cos 9) = e - aN ^ + ^ SL " (x)dx , (24) 

where N = J^°° n(x)dx is the total number of particles. Note that N is so large 
that exp (—(tN) — ► 0. Let us check different cases. We have 

cosfl>0 =*> P(-oo,cos0) -> 0, P(oo,cos0) = 1, 

(25) 

cos6»<0 => 7>(-oo,cos0) = 1 , PCoOjCos^) -> 0, 

The result for cosf? < 0, x = — oo seems to be a little confusing, but it just shows 
that a particle going in negative direction and starting from x = — oo has no particle 
to collide with. We are interested in the region with positive x: 

V(O,cos0) = e~ N °^ !1 + N '^ , (26) 

where N' — f_ n(x)dx. Let us assume that N and N' are such that 
exp (-aN') -> and exp (-ct(JV - AT')) -> 0. In this case 

7>(O,cos0) ->0, Vcos6>. (27) 

This condition is easily satisfied , because there is always a point A such that 
J^ 4 n(x)c?x = y, and then we should shift origo of our frame to this point along 
the x axis. 

This escape probability determines the free particle distribution ff ree (x,p) as 
a fraction of the total particle distribution 

ffree(x, P )=Vf(x,p), (28) 

and the interacting part of the particle distribution fi„ t is defined as 

f int (x,p) = (l-V)f(x,p), (29) 



Kinetic freeze out 



15 



where f(x,p) = ff ree {x,p) + fint( x >P) ■ The total density n(a;) is given as 

N°(x) = n(a;) = / d 3 pf(x,p). (30) 



Note that for x — > — oo, n(x) — > constant. Condition ( f27|) means that we do not 
have frozen out particles at x = 0, i.e., f(0,p) — fint{0,p). It is obvious from eq.(|24|) 
that if n is constant, then V becomes identically zero, because N — > oo, and the 
post-FO component can never emerge. Another important fact that could be seen 
from cq. (pi|) is that if 

7>(O,cos0) = 0, cos6»<0 

it will be equal for all x > 0. So backward going particles can not freeze out 
in our consideration. From this point of view this volume emission model close to 
idealized model with drain term discussed in section|]3. We will see later that the 
results of these calculations are similar in some aspects to those in the kinetic freeze 
out models discussed above. 

If the system is truly one-dimensional for all x values, then the total density 
n should vanish for large x, otherwise V vanishes. However, this contradicts the 
assumed conservation of flux in the stationary case: 

A^oc) = N^Q) = J d z p(p x /E)f(w,p) ? 

is incompatible with 

A°(oo) = n(oo) = J d 3 pf(oo,p) = 0, 

since the velocity p x /E < 1. Therefore, to get a stationary one-dimensional flow 
the system should have a finite size in the freeze out direction. 

Suppose that there exists a boundary at x = L > 0, so that for x > L, the 
density falls off very rapidly and the escape probability is almost zero there. Such a 
situation happens for a semi-infinite tube open to the vacuum at x — L. We should 
write 

N= I n{x)dx. (31) 

J — oo 

We are going to show that the equations ( p8| , |29| ) together with the conservation 
laws determine all the distribution functions when we assume a thermal spectrum 
for the interacting component. 

First, note that all distributions are specified if, 

1. the interacting flow velocity, Vi n t(x), 

2. the interacting temperature, T(x), 

3. the interacting density, ni„t(x), and 
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4. the escape probability in the x direction, 

P (x) = efo an{x)dx , (32) 
are known. To see this, first we write 

V{x,cos9) = P 1 {cose){P Q (x)} 7 ^ s , (33) 

where 

P^cosO) =V(0, cos 6) = e-^^^+^'s^ , (34) 
and express the total and free distributions in terms of fmt- 

f{x,p) = T ^fint(x,p) = T ^n i nt(x)^e- p "< t/T , (35) 

ffree(x,p) = T ^n i „ t ( a; )ie-^< t /T ) (36) 

where 

with 7 = 1/ yl — Vi n t{x) 2 and Z is the normalization factor, 

z = z{t) = J d 3 P e -p*«r/r 

In our stationary regime, the conservation laws are expressed as 

N 1 (x) = const — A fl (0), T m (x) = const = T 01 (0), T n (x) = const = T n (0), 

(37) 

where 

AT 1 ^) = J d 3 p^f(x,p), 
T 01 (x) = J d*pp x f{x lP ), 

T"(x)^ J d 3 p { -^ff(x,p). 

Once the initial values A^ 1 (0), T ol (0) and T 11 (0) are specified, these equations, 
together with Eq. (j35f), determine algebraically T(x), Vi n t{x) and rij n t(a;), at each 
x as functions of Pq(x). 

On the other hand, from Eq.([32"|) 
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and 

N°(x) = n(x) = J d 3 pf(x,p) = n mt (x) J ^-l_V^r/T 
so that we get an integro-differential equation for P , 

P dx an ^ x >] dp l- Pl ( C o S 6)P V™e Z € 
To compute n — N a , N 1 , T 01 and T 11 , we need to know the integrals: 
h[Pint,T,v int \ s |^ I -_^^ e^<VT | 

r m m i f ,q PCOS0 1 
^2 [Pint, T, Vint }= d*p — — — — 

7 po 1 - PxicosBjPo 1 / 



(39) 



(40) 



. e -P M <*/T 

i/ cos e e ' 



h[Pint,T,V int ]= j d3 P 1 _p i 

/' 
d Z V V - 



(cos^PoVcose 6 



• 2 cos 2 9 I - P "«r /r 



po 1 - Pi(cos6»)P 1/cosee 

Note that these functions arc not scalar, and the above expressions are valid in 
the frame where the density distributions are at rest, i.e. in the LR frame. The 
local rest frame quantities are labeled by *. Then in the local rest frame p — > p* , 
p^u^^p* = E*, and: 



l/cos0 



e -E*/T 



h [Pint , T, V mt ] = [ d 3 p* / 

J l-P 1 (cos9){P (x)} 1 

I 2 [P in t,T,v in t]= [d*p* 

J & 1 - Pi (cos ffj {Po(x)} ' 

I 3 [Pi nt ,T,v in t}= [d 3 p* ; P * m * 6 ] ^ 1/cose e- E ' /T 

J 1- P 1 {cos9){P (x)} 1/cos9 

Mp rp !_ f ,3 . P*Wfl* 1 _ 

I 4 [Pint,T,v mt ] = dp — — ————^e 

J & 1 - Pi (cost') {Po(X)\ ' 

with 

cost* =p x /p 



p* COS 6* 1 -E*/T 



P* 2 — p* 2 cos 2 9* + (jp* COS 9* + JV int E* 
and the limits of integral are restricted by jp* cosd* + "fVi nt E* > 0. 
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Fig. 5. The freeze out diatribut 
(RFF) calculated according to eqs 



ff ree (x,p), in the Rest Frame of the freeze out Front 
A, B and C correspond to x = 30fm, 80fm, 89.66fm 
respectively and u^ FG \ x=Q = (1,0, ~0j~ We put N = 180 , N' = 90 , m = , n = 1 and a = 

fm 2 = 10/ibarn. Contours correspond to ff ree = [1, 2, 4, 6]-10~ 35 (A), 10" 13 (B), 10" 9 (C*) ^— - 




McV 3 



Note that the form of the final freeze out distribution (C) seems to be similar to the cut Jiittner 
distribution for large p x . 



5.1. Extrapolation approach 

It is obvious that solving integro-differential equation (fio|), and evaluating the above 
integrals require nontrivial numerical calculations. Nevertheless, we would like to 
show a simple approach of the solution of this model. We can not get the equa- 
tions for 

^intipz) ,11^ {x) and T(x), since we can not evaluate the integrals in the 
conservation laws, but we are going to extrapolate ff ree (x) , f mt (x) and n(a;). 
Let us define a set of points on the interval [0, L] so that: 

{xi},i = 1, ...,m, Xi c [0, L] , xi — 0, x m = L . (41) 

Ax ■ 

A Xl = x l+ i - Xi, < < 1 . (42) 

We assume that / x (p) = f(0,p) = f{ nt {p) = f mt (0,p) is known (i.e. n(0),^'(0) 
and T(0) are known). The normalization factor Z can be evaluated as 

Z{T) = f d 3 pe- p "< t/T = 'lTTjT 3 a(2K 1 (a)+K (a))^l8TTjT 3 . (43) 

RFF J 

We have calculated Z(T) in RFF using the same mathematics as in the previ- 
ous section, and all calculations below will be made in RFF too. The probability 
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Pi (cos 8) = P(xi, cos 9) is then 

_ A , a e^g) +jv ,_^_ + _^_ g (n i+ n i+1 ) 

-P 4 (cos6») = e J=1 , Pi (cos 0) = 0. (44) 

If n(x) is a slowly varying function of x, we can extrapolate the escape probability 

by 

P m (cos0) - ^(cos^es***"*"* . (45) 

Since P, -C 1 for almost all values of Xj, except for the few last nodes, /*"' is a 
much more smooth function of x than ff ree i and its extrapolation will be: 

fU\ip)= (l-A+i(cos0))/,(p). (46) 



Next we calculate: 



/ i+ i(p)-/m(p) + // ree (p)- W 
fii+i - / d 3 p/m(p), (48) 



P m (cos0) = ^ (cos 6>) e^ Ax >^^- , (49) 
//+T(P) =^<+i(co B e)/ i+ i( P ), (50) 

/i + i(p) = /i?i(p)+//+T(p), (si) 

n,+i = / d 3 pf. l+1 (p). (52) 

Results received according to such an extrapolation are shown in the Figures 
| - 0. We can see that n(x) decreases very sharply and vanishes at point x = L, 
determined from the condition N — N' + n(x)dx. We can also observe the 
similarity of Figures || with the cut Juttner distribution. At the same time the 
initial stages are much more elongated in the x direction in Figures || and [| what 
may also lead to a large-pt enhancement. This model leads to incomplete freeze out 
as the one, presented in sectiorfj}., but such a problem can be cured if rescattering 
would also be included in the model, as we demonstrated in the previous case. 

So far in the calculations described above we did not take into account the 
conservation laws. Let us check do we break conservation laws or not. According 
to (||, ||, |l|) we have: 
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Fig. 6. The same as Fig. g just for w^ FG |x=0 = 7(1, 0.5, 0, 0). 

/i+i(p) = (l-Pi +1 (cosfl))/ i (p)+P i+1 (cosd) [(l-P i+ i(cos0)) fi(p)+ 

+P 4 [(l -P 4 (cos0)) fi-x(p)+Vi-i [(l - 7Vi(cos0)) /^(p) + 

/<+i(p) ~ /<(P) = (^+i(cos0) - P i+1 (cosd)) + O (V 2 ) , 
since P 4 < 1 as well as P; < 1. Using (||, ft9|) we get: 

/i+i(p) - /<(P) = P 4 (cos0)e^ A ^ ( e ^ A ^ - 
Let us define 

jr. s rf3 P x d3 P x x 

dF = —f-y or — —p x p u or —f-p p . 

pV pV pU 

Then the conservation laws take the form 



J dF fi(p) = Const, 



(53) 



(54) 



/i(p). (55) 



(56) 



where Const = A fl (0) , or T ol (0) , or T n (0) respectively. Assume that this is true 
for fi{p), and then let us check 



Vi(cos6) = Pi (cos 0) 



(57) 
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Fig. 7. Total density n(x).We assume N = 180, N' = 90, m = , no = 1 and a = fm 2 = 
10 /ibarn. We can see n(x) is almost the same for x C [0, 87] fm, but then it decreases very sharply. 



where n.; = ^ ^- x i "' + , and Pi (cos 6) is defined by ( p4| ) . The changes of 
j'=i 

conserved values at the ith step of extrapolation are given by 

dFif i+l{p )-u P ))=[ U + ^Y^^^Hv)^e)e^ d F. 

J \ cos ft 2 / cos ft 2 

(58) 

Since nj +1 rs n i+1 with a good accuracy, and Vi(cos9) has its maximum at cos# = 
±1 (if we put N' = f we get 7^(1) = Pi(-l)) we have 

| dF(f i+1 (p) - ^(p)) < Pi(l) e^ a - (Ax,) 2 • Const. (59) 

Finally, 

J dF (f k (p) - Mp)) < Px(l) Const °- YjF* (Axi) 2 . (60) 

We have seen that ( |£) . < 1 (Fig. 0) and 7>i(l) < 1 for all i, except the last few 
nodes. So, making Ax^ small enough, we can keep dN x , dT 0x , dT xx conserved 
with the necessary accuracy. 
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6. Conclusions 

In this work we evaluated in a simple kinetic model the freeze out distribution, 
ffree(x,p), for stationary freeze out across a surface with spacelike normal vector, 
da^da^ < 0. 

The first simple kinetic freeze out model, adopted from [|| (see section |3.) 
reproduces the cut Juttner distribution as the limiting distribution, // re e, after 
complete freeze out at large distances. However, the model at the same time leads 
to unrealistic consequences, namely that the interacting part of the distribution, 
fi n t also survives fully, as the other part of the Juttner distribution. Thus, having 
both components at the end in this model, the physical freeze out is actually not 
realized. 

Here we have presented a solution for an improved but still rather approximate 
kinetic freeze out model which takes rescattering into account (see sectio|]4.). In this 
model the interacting component is assumed to be instantly re-thermalized taking 
a spherical Juttner shape at each time step with changing parameters. The three 
parameters of the interacting component, fi n t, are obtained in each time step. The 
density of the interacting component gradually decreases and disappears, the flow 
velocity also decreases and the energy density decreases also. The temperature, as a 
consequence of the gradual change in the emission mechanism, gradually decreases 
at the final stages of the freeze out, when only high energy, forward going particles 
are taken away from the interacting component. 

The arising post freeze out distribution, ff ree is a superposition of cut Juttner 
type of components, from a series of gradually slowing down Juttner distributions. 
This leads to a final momentum distribution with a more dominant peak and a 
forward halo, Fig. ||. In this rough model a large fraction (~ 95%) of the matter is 
frozen out by x = 3A, thus the distribution // ree at this distance can be considered 
as a first estimation of the post freeze out distribution. One should also keep in 
mind that the model presented here does not have realistic behavior in the limit 
x — > oo, due to its one dimensional character. Nevertheless, this improved model 
with rescattering enables complete freeze out. 

The second model presented here (see sectioi£]5.) shows that we can achieve full 
freeze out in finite length even in oversimplified one dimensional models. Thus, the 
drawback of the previous model that full freeze out happened exponentially slowly, 
can be remedied by including more realistic features in the model in a straightfor- 
ward way. 

These studies indicate that more attention should be paid to the final freeze 
out process, because a more realistic freeze out description may lead to large pt 
enhancement [|i],|25| as the considerations above indicate (Fig. |^). If the heavy 
ion reaction is basically described by a kinetic model of weakly interacting hadrons, 
then idealized FO models with an assumed FO hypersurface are on the limits of 
their applicability, and even then, such models could only be applicable for the 
heaviest systems. 

If, however, QGP is formed in heavy ion reactions, the number of degrees of 
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freedom increases tremendously, and continuum models become more suitable to 
the problem than dilute gas kinetic or string models with binary interactions. In 
case of rapid hadronization of QGP and simultaneous freeze out, the idealization 
of a freeze out hypersurface may be justified, however, an accurate determination 
of the post freeze out hadron momentum distribution would require a nontrivial 
dynamical calculation. 
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Notes 

a. For the classification scheme see jittp : / / www . aip . org/pubscr vs/ pacs . htm] . 

b. Assuming that the matter is characterized with 4- velocity Uj, f q, which is nor- 
mal to the freeze out surface in the three dimensional space, and differs from 
the Local Rest frame (LR) velocities of f FO , (i.e., vf^ and u F ) we introduce here 
n(p,T) = 87rT 3 e^ /T (27r?i)- 3 , a=f, so that n{pb,T) = na 2 K 2 (a)/2 is the in- 
variant scalar density of the symmetric massless Jiittner gas, b = a/y/1 — v 2 = 
07, v = = dao/dcr x = u rfg \rff, A = (2 + 26 + b 2 )e~ b , and 

Kn(z,w) = z- n Tdx (x 2 - z 2 )^ 1 ' 2 e~* , 

(Zn). J w 

i.e. K n (z,z) — K n (z). When evaluating the limit we used the relation 

JC n (a,b) a ^—? K n (a) (n — l)!2 rl_1 a _ ™. This baryon current may then 
be Lorentz transformed into the Eckart Local Rest (ELR) frame of the post FO 
matter, which moves with u% = N^/{N v N v f/ 2 = 1e{1,v Ei 0,0)\rfg in the 
RFC, or alternatively into the Rest Frame of the Freeze out front (RFF), where 
da^ = (0, 1, 0, 0)\rff and the velocity of the RFC is u RFG = 7^(1, v, 0, 0)\rff- 
Then the Eckart flow velocity of the matter represented by the cut Jiittner 
distribution viewed from the RFF is u F = 7 C (1, v c , 0, 0)\rfFi where v c = 
(v + v E )/(l + vv E )- 

c. We are considering fast particles, so that v re i ~ v. A formula similar to ( |24| ) is 
also obtained for massless particles. 
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